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Abstract 

We study future observational constraints on cosmic string parameters from various types of 
next-generation experiments: direct detection of gravitational waves (GWs), pulsar timing ar- 
ray, and the cosmic microwave background (CMB). We consider both GW burst and stochastic 
GW background searches by ground- and space-based interferometers as well as GW back- 
ground detection in pulsar timing experiments. We also consider cosmic string contributions to 
the CMB temperature and polarization anisotropies. These different types of observations offer 
independent probes of cosmic strings and may enable us to investigate cosmic string properties 
if the signature is detected. In this paper, we evaluate the power of future experiments to con- 
strain cosmic string parameters, such as the string tension G\i, the initial loop size a, and the 
reconnection probability p, by performing Fisher information matrix calculations. We find that 
combining the information from the different types of observations breaks parameter degenera- 
cies and provides more stringent constraints on the parameters. We also find future space-borne 
interferometers independently provide a highly precise determination of the parameters. 



1 Introduction 



Cosmic strings are linear topological defects which is formed at spontaneous symmetry breaking 
(SSB) in the early Universe pQ (as a review, see Ref. [2]). Some inflation models based on super- 
string theory predict fundamental strings, D-strings and their bound states of cosmological length. 
Such stringy cosmic strings are called as cosmic superstrings [3H5]. They form a complicated string 
network, which consists of infinite strings and closed loops, and may leave remarkable signatures in 
the universe through their nonlinear evolution. If their signals are observed, not only the existence 
of cosmic strings will be confirmed but their properties might be studied. This enables us to ob- 
tain implications on both the history of the universe such as inflation or SSB, and physics beyond 
the Standard Model of particle physics such as grand unified theory (GUT) or superstring theory 
through the study of cosmic strings. Therefore, cosmic strings are important probes on both cosmol- 
ogy and particle physics, and we are motivated to study how the properties of cosmic strings can be 
determined by future experiments. 

Various types of observational signatures of comic strings have been studied intensively. One of 
them is the gravitational wave (GW) [BlTT?]. The main source of GWs in the string network is cusps 
on loops. A cusp is a highly Lorentz boosted region on a loop which appears 0(1) times in an 
oscillation period of the loop and it emits a strong beam of GWs, which we call a "GW burst". 
The GW bursts can be detected in two different forms. One is "rare bursts", which are infrequent 
but strong enough to be detected alone. The other is a stochastic GW background, which consists 
of many small bursts overlapping each other [TUlfTT]. In our previous paper [2D], we studied how 
cosmic string parameters are determined by direct detection of GWs in future ground-based GW 
experiments, such as Advanced LIGO [21], Advanced Virgo [22] and KAGRA [23], and showed that 
measurements of the burst rate and the GW background provide different information and leads to 
better constraints on parameters when they are combined. 

In this paper, we extend our previous study in two ways. First, in addition to the ground- 
based detectors, we consider other types of GW experiments, that is, direct detection by space-borne 
interferometers such as eLISA/NGO [21], BBO [25] or DECIGO [26], and observation of the GW 
background in pulsar timing experiments, such as Parkes PTA [27], EPTA [28], NANOGrav [291130] 
or SKA [31]. These different types of GW experiments provide different information, since each type 
of experiments has its best sensitivity at different frequency: ~ 10 2 Hz for ground-based detectors, 
~ 10~ 2 Hz for eLISA/NGO, ~ lO^Hz for DECIGO and BBO, and ~ 10~ 8 Hz for pulsar timing 
experiments. Gravitational waves in different frequency bands are emitted at different redshifts and 
carry information of cosmic strings living in different epochs of the Universe. 

Second, we also take into account observation of the cosmic microwave background (CMB). 
Cosmic strings may induce temperature and polarization fluctuations in CMB through gravitational 
effects, which is also an important observational signature of cosmic strings and studied intensively 
(for example, see Refs. [3"2THo] ). Although current CMB observations indicate that cosmic strings are 
not the dominant source of CMB fluctuations, the string signatures might be observed in small scale 
fluctuations and/or B-mode polarization by future experiments such as Planck [UJ and CMBpol 

1 GWs from kinks may dominate that from cusps in the case of cosmic superstring [T5J|TS] . However their con- 
tribution strongly depends on the fraction of loops with junctions. So we do not consider their contribution in this 
paper. 
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|48j . Expected constraints on cosmic string parameters by these experiments are already studied in 
Ref. [62J by Fisher information matrix calculations. In this paper, we consider combination of the 
constraints with GW experiments, which are expected to provide different information on cosmic 
strings. 

As in Ref. [20], in this paper, we focus on three parameters which characterize cosmic string 
network. The first one is string tension fi, or the product of it and Newton constant G, G\x. It 
represents the energy stored per unit length in a cosmic string. For the field theoretic string, it is 
comparable to the square of the energy scale of SSB, while for the cosmic superstring, it is determined 
by the energy scale of the superstring theory and the warp factor of the extra dimension where the 
string is located. The value of Gfi affects not only the amplitudes of the GW bursts and GW 
background but also the spectral shapes through the change of the lifetime of loops. The amplitude 
of CMB fluctuations is also affected by the value of G[i. 

The second one is initial loop size a. The typical size of a loop at its formation is characterized 
by at, where t is the time of the loop formation. In principle, the value of a can be predicted, if we 
can solve the nonlinear evolution of the string network, and there are many works which attempt to 
determine the value of a with numerical or analytical methods [19H60"] . However, it is not yet clearly 
understood, so we treat a as a free parameter. The value of a affects both the GW burst rate and 
the GW background spectrum. On the other hand, the CMB signature is independent from a, since 
it is induced mainly by infinite strings and the contribution from loops is negligible. 

The third one is reconnection probability p. For field theoretic strings, p is roughly equal to unity 
while for cosmic superstrings, it can be much smaller than 1. The string network becomes denser 
as p gets smaller. This leads to the enhancement of the burst rate, the GW background spectrum, 
and the CMB fluctuations. The shape of the CMB power spectrum is also affected, since the typical 
length scale of the string network becomes smaller and the average velocity of the strings become 
larger. 

In this paper, we aim to investigate the abilities of the different types of experiments to constrain 
string parameters and study how they complement each other. We calculate the burst rate and 
the GW background spectrum using the formulation described in Ref. [2D]. For CMB, we use 
CMB ACT [61], which calculate the power spectrum by approximating the string network as an 
ensemble of randomly oriented straight segments [35] . Finally, we predict constraints on Gfi, a and 
p from the future experiments by performing Fisher matrix calculations and find that the different 
types of experiments break the degeneracies in the parameters and help to tighten the constraints 
when they are combined. 

This paper is constructed as follows. In Sec. 2, we briefly review the model to describe the cosmic 
string network and the methods to calculate GWs from string loops and CMB power spectra induced 
by strings. In Sec. 3, we describe the sensitivities of current and future experiments and the Fisher 
matrix formalism for direct detection, pulsar timing, and CMB experiments. In Sec. 4, we first 
show the parameter space to be explored by these experiments, and then calculate constraints on the 
cosmic string parameters for three fiducial models. We summarize the paper in Sec. 5. Throughout 
the paper, we use the cosmological parameters from the 7-year WMAP data (WMAP+BAO+H0 
mean) [67]: the ratio of the present energy density of baryon to the critical density f2&/i 2 = 0.02255, 
that of cold dark matter Q c h 2 = 0.1126, that of dark energy Q\ = 0.725, the spectral index of 
the primordial curvature perturbation n s = 0.968, the reionization optical depth r = 0.088, the 
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amplitude of the primordial curvature perturbation A^(fco) = 2.43 x 10 9 at k = 0.002Mpc 
Hubble constant H Q = 70.2km/s/Mpc and the primordial helium abundance Y p = 0.326. 



2 Calculation of various types of cosmic string signatures 

In this section, we briefly mention how to calculate the observational signatures of cosmic strings. 
After explaining the analytic model of the cosmic string network, we describe the formalism to 
calculate the burst rate, the GW background spectrum and the CMB power spectrum and mention 
their dependence on string parameters. 



2.1 The model of the cosmic string network 



As in our previous paper [20], we adopt the model in Refs. , which is based on the velocity- 

dependent one-scale model [65J. The network of infinite strings is considered as a random walk with 
a correlation length £, which corresponds to the typical curvature radius and interval of infinite 
strings, and the total length L of infinite strings in volume V is given by L = U/£ 2 . The equations 
for 7 = £/t and the root mean square velocity of infinite strings v are given by 

^ = -i + m + ^ + mv\ (i) 

7 at 27 
dV -'l- v 2 )H ( k M - 2v) , (2) 



dt \Ht-y 
where k(v) = jnif^ [66] - The Hubble parameter is given by 

H(t) = - = H [n A + fi m (l + zf + Q r (l + zf] 1/2 , (3) 
a 

where a(t) is the scale factor, 1 + z — ao/a(t) is the redshift, ao is the present value of the scale 
factor, Q m = Qb + and Q r is the ratio of the present energy density of radiation to the critical 
density. The parameter c(t) represents the efficiency of loop formation. The value of this parameter 
between the radiation-dominated era and the matter-dominated era is interpolated by 

m = TT ^f, (4) 

where we set c r = 0.23, c m = 0.18 and g = 300 [351165] . 

The parameters 7(2) and v(t) have different time evolutions depending on the Hubble expansion 
rate. In our previous paper, we used the asymptotic values for each radiation-dominated and matter- 
dominated Universe. However, in this paper, since we additionally investigate the effect on the CMB 
fluctuations, which are produced near the transition from the radiation-dominated to the matter- 
dominated Universe, we numerically solve Eqs. ([T]) and ([2]) to evaluate the values of 7(2) and v(t). 
Their initial values are determined by the solution of d'y/dt = and dv/dt = in the radiation- 
dominated era with c = c r and Ht = 1/2. 
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(a) 7 for p = 1, 1CT\ 1CT 2 . (b) v for p = 1, 1CT 1 , 1CT 2 . 

Figure 1: The time evolution of 7 and v for different values of p. 



Figs. 1(a) and 1(b) show the time evolutions of 7(2) and v(t) and their dependence on the 
reconnection probability p. For p = 1, 7 and t> reach the asymptotic values in the matter-dominated 
era from those in the radiation-dominated era by z ~ 10 2 . After the dark energy becomes dominant 
component of the universe at z < 1, 7 begins to increase and v begins to decrease, because the 
exponential expansion of the universe dilutes strings and makes them slow down. For small p, we 
see the overall magnitude of 7 decreases as p decreases, the dependence is proportional to p' 1 and 
p" 1 / 2 for the radiation- and matter-dominated era, respectively [20J. The asymptotic value of v is 
approximately 1/V2 in both the radiation- and matter-dominated eras. 



2.2 Gravitational waves from cosmic string loops 

Here, we briefly describe the formalism to calculate GWs from cosmic strings. The detail is described 
in our previous paper [20], or originally in Refs. p7QHT3] . 

Since the main source of GWs in the string network is cusps on loops, we first evaluate the 
density of loops in the Universe in order to calculate GWs from them. In the scaling regime, where 
the typical length scale of infinite strings is proportional to the Hubble scale and their number in a 
Hubble horizon remains constant, infinite strings continuously convert their length into loops. The 
number density of loops formed at time ti is given by 

dn (++\j+ dti { aifi)\ 3 

dt^ ti)dti - ary(U)*t}\a(t)J ' (5j 

at time t. Loops continue to shrink by releasing energy as GWs and eventually evaporate. The 
length of a loop formed at t{ is given by 

l(t,U) =aU-TGfx(t-U), (6) 

for time t. Here, T is a constant which represents the efficiency of the GW emission from loops and 
set to 50 in this paper. If a 3> TG[i, loops are long-lived, that is, they survive more than a Hubble 
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time. On the other hand, if a <C TGfi, loops are short-lived, and they evaporate within a Hubble 
time. Using Eqs. (j5J) and (EJ), we can describe the number density of loops in terms of length I and 
time t. 

Cusps emit GWs of / 3> Z" 1 into a small solid angle. The linearly polarized waveform of a GW 
burst emitted in a direction n by a loop with length / at redshift z is given by 

h^t, n) = J dfh(f, z, /)e- 2 ^ t e+ l/ (n) x 6(n ■ n c - cos [9 m (f, z, /)]) x 6(1 - 9 m (f, z, I)), (7) 

where n c is the direction of the center of the burst, 9 m is the beaming angle of the GW burst which 
is given by 

6 m (f,zJ) = ((l + z)fl)- 1/ \ (8) 

and e+ = l fl m u — l u m fl is the polarization tensor (for plus polarization), where 1^ = (0, 1), = (0, m) 
and 1 and m are unit vectors orthogonal to n and each other. The Fourier transform of the GW 
amplitude, h(f,z,l), is given by 

H f, zJ) ^2M w -^L m] , (9) 

where r(z) = J* dz' /H(z'). The first Heviside step function in Eq. ([7]) is introduced to account 
for the beaming of the GW and the second one is for the low frequency cutoff at / < 

The arrival rate of GW bursts with frequency / and amplitude h emitted at redshift z is given 

by 

1 (f,h,z) = ty m tf,z,i) ,°, u * J4 * (^i) 3d ^e(i-e m (f,zM (io) 



dzdh KJ ' 11 4 mw ' ' J (1 + z^-fit^atja + TGfi \a(t) J dz 
where 



dV 4^{zy{z) 
dz [ ] ~ H(z)(l + z)- [ } 

From Eqs. ([6]) and (Q, I and tj can be expressed as 

a + 1 G/i 

which enables to express Eq. ( Tl0|) in terms of h, z, and /. Finally, the total arrival rate of GWs 
today for given frequency and amplitude is given by 

dR _ r°° dR 
dh J dhdz 

The GW bursts are identified as a single burst if they do not overlap each other. In contrast, 
bursts overlapping each other are observed as a GW background. Namely, a GW background is 
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formed by bursts which come to the observer with a time interval shorter than the oscillation period 
of themselves. According to the criteria in Ref. [13] . such bursts have amplitude smaller than h*, 
which is determined for a given frequency as 



Then, the amplitude of a GW background, ^gw(/) = (dpcw/d In f)/p cr , where pew is the energy 
density of the GWs and p cr is the critical density of the universe, is given by 



The main contribution to Qqw at frequency / comes from loops expiring at the redshift which 
satisfies min{a, TGfi} xt~ + z)' 1 in the radiation-dominated era or loops expiring recently 

Here, we briefly mention the parametric dependence of the burst rate and flow- The burst rate 
and the amplitude of the GW background are basically enhanced as G\x increases. However, the 
value of Gfi also affects the spectral shape of the burst rate and GW background, since the number 
density of loops are affected by G/i through the lifetime of loops. This is important especially in 
the case of Gfi > a, where G/x determines the size of expiring loops. The large value of a decreases 
the initial number density of loops. However, it also has an effect to increase the number density of 
loops, since large a makes the lifetime of loops longer. Therefore, the effects of a on the burst rate 
and the background depend on values of /, h and the other parameters. Larger p simply leads to a 
larger burst rate and larger amplitude of the background, through the factor 7" 2 in Eq. (0E). We 
refer to Ref. [20] for more detailed discussion. 

2.3 CMB fluctuations induced by cosmic strings 

Here, we describe the method to calculate the power spectra of CMB fluctuations produced by cosmic 
strings. We refer to Refs. J33EI] for the details. 

The evolution of the string network is highly non-linear and produce all types of perturbations: 
scalar, vector and tensor modes. For the calculation of the CMB power spectra, we use a code based 
on CMBACT [61], which we modify to include the time evolution of 7 and v. This code is based on 
the semi-analytical method described in Ref. [35J, which models the string network as an ensemble 
of discrete straight line segments. The length and the velocity of each segment are set to the solution 
of Eqs. (JTJ) and (j2J). The position and the direction of the velocity are randomly selected for each 
segment. At each time step, some segments are removed so that the number density of strings is 
consistent with the scaling. Then, we can derive the energy-momentum tensor of such a simplified 
network and compute the power spectra of CMB fluctuations. In this paper, we assume that infinite 
strings have no wiggliness, which is introduced in Ref. [35]). 

Here, we briefly discuss the parametric dependencies of the CMB power spectra. More detailed 
discussion is given in Ref. [H] . In Figs. [2] and [31 we show the CMB power spectra of temperature 
fluctuations cf T,str and B-mode polarization cf B ' str induced by cosmic strings, respectively, for var- 
ious parameter sets. We also show the temperature spectrum induced by the inflationary primordial 
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(a) Dependence of the temperature spectrum on (b) Dependence of the temperature spectrum on 

Gfi. The green, black and red solid lines represent p. The black, red and green solid lines represent 

the spectra for Gn = 10" 6 ' 5 , Gfx = 1(T 7 and Gfi = the spectra for p = 1, p = lCT 05 and p = 10" 1 , 

10~ 7 - 5 , respectively, with p fixed to 1. respectively, with Gfi fixed to 1CT 7 . 



Figure 2: The power spectra of CMB temperature fluctuations induced by strings for various param- 
eter sets. In each figure, the black solid line represents the spectrum for the fiducial parameter set 
taken for the Fisher analysis in Sec. 4.3. We also show the spectrum predicted by inflation as a blue 
dotted line and the expected noise levels for Planck and CMBpol as a black and orange dotted line, 
respectively. 



perturbations cJ T,m * in Fig. [2j the B-mode polarization from the inflationary GWs C BB ' m ^ and that 
from the gravitational lensing of E-modes c BB,len in Fig. [3], which are calculated by CAMB [921193] . 
In addition, we plot the expected noise levels of Planck and CMBpol, which is given by the sum of 
the instrumental noise and the cosmic variance 



fc^iV™ = !£±iU -A- -{ C T™ S + N T ,), (17) 
2tt t,/ 2tt y (2/+ l)/ v 1 ' h v ' 

for temperature fluctuations, and 

l l±Q.N^i = iii±Hj-A_(Cf^ en + N P ,), (18) 
2vr b ' 1 2tt y (2/ + l)/ v 1 F ' lh K 1 

for B-modes. The instrumental noise spectra for the temperature Nt-,i and the polarization Np-i 
are defined in Sec. 3. 

The power spectrum of CMB temperature fluctuations induced by strings has a single bump, 
unlike the acoustic oscillations seen in the inflationary CMB power spectra. This is because strings 



The variance of CI is given by (AC?) 2 = ^(Cj + N a ,if, where i denotes TT or BB and a denotes T or P. 
When we take a logarithmically homogeneous binning of / with bin width A hi I = 1, there are I multipoles in a bin at 
I. Since different multipoles arc independent, the noise level par each bin should be given by ACf/Vl- This allows us 
rough estimation of the detectability of the signal by eye in Figs. [2] and [3] 
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(a) Dependence of the B-mode spectrum on G\x. (b) Dependence of the B-mode spectrum on p. The 
The green, black and red solid lines represent the black, red and green solid lines represent the spec- 
spectra for Gfi = 10~ 6,5 , G\i = 10 -7 and G\i = tra for p = 1, p = f0~ 05 and p = 10 _1 , respec- 
10~ 7 5 , respectively, with p fixed to 1. tively, with G\i fixed to 10~ 7 . 



Figure 3: The power spectra of CMB B-mode polarization induced by strings for various parameter 
sets. In each figure, the black solid line represents the spectrum for the fiducial parameter set taken 
in the Fisher analysis in Sec. 4.3. The blue dotted line is the spectrum induced by the inflationary 
GWs. The tensor-to-scalar ratio r is set to 0.1. The purple dotted line is the spectrum of B-mode 
polarization generated through the gravitational lensing of the inflationary E-mode polarization. The 
black and orange dotted line represent the expected noise levels for Planck and CMBpol, respectively. 



generate fluctuations constantly and such fluctuations are incoherent, while the primordial fluctu- 
ations generated during inflation oscillate coherently. The peak of cJ T,str corresponds to the scale 
of the perturbations generated at last scattering. Another important difference is that cJ T)Str does 
not decay exponentially as I increases in contrast to the Silk dumping of the inflationary spectrum 
at large I. This makes cJ T ' str larger than cf T,m ^ at high / and it may be observable by future 
experiments^]. The B-mode power spectrum also does not have oscillations, but has two bumps. The 
bump at high I and that at low / correspond to the polarization generated around last scattering 
and around reionization, respectively. If the amplitude of the primordial gravitational wave is small, 
strings can be the main source of B-mode polarization. 

The dependence of CMB power spectra on Gfi is simple. The amplitudes of temperature and 
polarization fluctuations are proportional to Gfi, so the power spectra are proportional to (Gfi) 2 . 
The effect of p arises through the correlation length of the string network 7 and the r.m.s. velocity 
of strings v. For small p, 7 becomes smaller and v becomes larger. If p is small, the small value of 
7 makes the string network denser and enhances the amplitude of the CMB fluctuations. Also, the 
small correlation length 7 makes the typical scale of perturbations smaller. However, at the same 
time, large v makes the typical scale of perturbations larger j5T]. Since the two effects on the typical 
scale of perturbations compensate, we find no apparent shift of the peak in Figs. [2] and [3] with the 



3 The power spectra of E-mode polarization C EE and the cross-correlation between temperature and E-mode Cj E 
by strings also decay more slowly than the inflationary spectrum at high I. However, they have small amplitude and 
do not affect the observation when we assume realistic parameters. 
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change of the value of p. While these effects tend to compensate, small effects are seen in the width 
of the peak and the slope of the small-scale spectrum. 



3 Formalism to evaluate the sensitivity and the constraints 
on parameters in each experiment 

In this section, we describe the method to calculate the Fisher information matrix, which is used 
to estimate the power of each experiment to determine the string parameters. For GW burst and 
background detection by interferometers, Ref. [20] provides more detailed derivations. 

Given a data-set of experiments, the parameter values that are most likely to result in the model 
prediction are those which maximize the likelihood function C. Such a method for parameter es- 
timation is called the maximum likelihood method and widely used in the analysis of cosmological 
observations [70] . The error in this estimation can be predicted by calculating the Fisher information 
matrix 

d 2 ln£ 

7lm - ~aw^' (19) 

Assuming a Gaussian likelihood, the expected error in the parameter Q\ is given by 



00, 



(20) 



3.1 GW burst detection 

The search for GW bursts signals from cosmic strings is performed by matched filtering [681169]. 
where we assume that the burst signal from a cosmic string cusp is linearly polarized and have the 
frequency dependence of /~ 4 / 3 . Then, the spectrum form is expressed by 

h+(f) = Af-^e(f h - f)Q(f - /,), (21) 

where the amplitude A can be read from Eq. (Q. The low frequency cutoff /; is given by the low 
frequency limit of the experiment and the high frequency cutoff fh is typically given by the most 
sensitive frequency of the detector [71] . The signal to noise ratio(SNR) p is given by 




where the GW signal h(f) = F + h + (f) is given by multiplying the detector response to plus polarized 
GWs F + . Effectively, F + can be replaced with the all sky- averaged value for orthogonal arm detec- 
tors, F+ ~ 1/V5for a single detector, or F + ~ 1 for the GW detector network which has 100% visibil- 
ity over the whole sky. The noise spectral density S n (f) is defined by (n(f)*n(f)) = S n (f)S(f — f')/2, 
where n(f) is the Fourier transform of the detector noise n(t) and (■ • • ) denotes the ensemble average. 
In this paper, we take the detection threshold as p > 4 [TO] . 
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For current LIGO, the noise spectrum is given by 

S n (f) = 1.09 x lO" 41 + 1-44 x lO" 45 



^30Hz\ 28 

1.28 x 10~ 46 ( 1 + ( ) ] Hz" 1 . (23) 



/90Hz\ 

w 




and we take / = 40Hz and fh = 150Hz [68J. The detection threshold p > 4 corresponds to the 
detection limit A ~ 9.1 x 10 _21 s -1 / 3 or fh ~ 1.7 x 10~ 21 at / = 150Hz. For the future ground-based 
GW detectors, such as Advanced-LIGO, we use 



Sn(f) = 10 



-49 



x _4.14 _ - + HI 



5 . /2-2x 2 + x 4 



x 2 V 2 + x 2 



Hz" 1 , (24) 



where x = //(215Hz) [72], and we take /; = 10Hz and fh = 220Hz. If we consider a world-wide 
detector network and assume F + ~ 1, the detection limit is A ~ 2.1 x 10- 22 s -1 / 3 or fh ~ 3.4 x 10 -23 
at / = 220Hz. For eLISA/NGO, we use 

Sn(f) = gO ^ + ^ + A™ ( 1 + ( f_ 




L 2 \ U41(£ 

where 5 acc = 2.13 x 10- (l + ^) m 2 S - 4 Hz- 1 , ^ sn = 5.25 x lO-n^Hz" 1 and ^ omn = 6.28 x 

10- 23 m 2 Hz _1 with the arm length L = 1.0 X 10 6 km j23], a take /, = 10~ 6 Hz and / h = 7.0 x 10~ 3 Hz. 
Then, the detection limit is A ~ 8.9 x 10^ 22 s^ 1/3 or //i ~ 5.8 x 10^ 21 at / = 7.0 x 10~ 3 Hz. For 
BBO or DECIGO, whose sensitivities are roughly the same, we use the configuration of BBO. It is 
designed to use a technique called time-delay interferometry, and the noise spectrum is given by 

Qm (RaSJ) , Mi) , Mf) Y 1 

s '' {f) ^{jM) + ^U) + ^U)J ' (26) 

The subscripts (A,E,T) denote the TDI variables and each noise spectrum is given by 

S A (f) = S E (f) = 8sin 2 (//2)[(2 + cos/> shot + 2(3 + 2 cos / + cos(2/)) 1 S acccl ], 

S T (f) = 2[l + 2cos/] 2 [S' shot +4sin 2 (/72)5 accel ], (27) 

where ^ shot = 2.0 x 10^ 40 /(L/km)- 2 Hz- 1 and S accc i = 9.0 x 10- 40 /(27r//Hz)- 4 /(2L/km)- 2 Hz- 1 with 
the arm length L = 5.0 x 10 4 km. For the calculation of the detector response Ra,e,t, see Ref. [73] • We 
take fi = 0.1Hz and fh = 0.25Hz, where the low frequency cutoff is determined to take into account 
the confusion noise from white dwarf binaries. This leads to the detection limit A ~ 1.2 X 10~ 24 s -1 / 3 
or fh ~ 1.9 x 10— at / = 0.25Hz. 

Let us suppose that we detect a sufficient number of GW bursts at /best, at which the detector 
is most sensitive. Under the assumption that the number of bursts follows a Poisson distribution 
[741476] , the Fisher matrix is given by 

f°° 3$ a$ i , 

where 3>(h) = dR/dh x T, which is a function of the model parameters, T is the observation time 
and /i m in is the smallest amplitude of the detectable burst at / = /best- 
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3.2 Search for the stochastic GW background by interferometers 



The GW background is searched by correlating output signals of two or multiple interferometers. 
The SNR in such a correlation analysis with N detectors is given by [771178] 



N N 



EE A 



1=1 j<i 



1/2 



where 



PU = - 5 V2T 



'If 



S h (f) 2 hu(f)\'' 
Ru(f) 



1/2 



(29) 



(30) 



Ru(f) = ( fj s h (f) 2 (| 7 /j(/)I 2 + mijj) + \s h {f) (nnU)SnAf) + ujSnAf)) + s n>I (f)s n ,j(f), 

(31) 

and T is the observation time. The subscripts I and J refer to independent signals obtained at 
each detector or TDI variables (A,E,T). The overlap reduction function between the 7-th and J-th 
detector ju is given by 



7lJ (f) = — / dn(F;(f,n)F^(f,n) + F^(f,n)F^(f,n))e 



-27ri/n-(xx-xj) 



(32) 



where Ff or Ff is the detector response to plus or cross polarized GWs of the 7-th detector, x/ 
is the position of the 7-th detector and Cl is the direction of GWs. We calculate 7/j following the 
procedure given in Ref. [79] for ground-based detectors and Ref. [80J for BBO. The signal spectrum 
can be corresponded to ficw as 

S h (f) = ^ficw(/). (33) 
In the weak signal approximation, Sh(f) S n j(f), Eq. (|30|) reduces to 



Pu 



3HI 
IOtt 2 



df 



|7/ J (/)| 2 ^gw(/) 2 
f 6 S nJ (f)S n M) 



1/2 



(34) 



The weak signal approximation is valid for p < 200 [78] [82]. 

From the non-detection of GW background, we can put an upper limit on Qqw- Current LIGO 
detectors set an upper bound of flew < 7.2 x 10~ 6 for 41.5Hz < f < 169.25Hz, assuming a flat 
spectrum. A 3-year run of the future detector network including Advanced LIGO, Advanced Virgo 
and KAGRA, would reach fi GW = 4.5 x 10~ 9 at 10Hz < f < 200Hz. A 3-year run of BBO/DECIGO 
would provide fi G w = 9.2 x 10~ 17 for 0.1Hz < f < 10Hz. Since eLISA/NGO is designed to have 
only one independent channel, we do not consider the cross correlation analysis in eLISA/NGO. 

The Fisher matrix for the GW background measurement is generally given by 



i m 



N N 

EE 

i=i j<i 



8T 
25 



df 



hu(f)\ 2 deMf)de m S h (f) 
Ru(f) 



(35) 
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Under the weak signal approximation, it reduces to 

= U?J 2T 2- I. X * rs^wwn ■ (36) 

3.3 Search for the stochastic GW background in pulsar timing experi- 
ments 

Pulsar timing experiments provide a unique opportunity to observe GWs in low-frequency band 
10~ 9 — 10~ 7 Hz [83H85] (for a review, see Ref. [86]). The analysis is based on the measurement of 
pulse time-of-arrival (TOA) variations. The stochastic GW background causes fluctuations in the 
TOAs. One can extract the signal from noise associated with individual pulsars by correlating TOAs 
between different pulsars. 

We follow the formalism described in Ref. [89] to calculate the SNR and the Fisher matrix for 
detection of GWs in pulsar timing experiments. Let us assume observations of M ^> 1 pulsars at 
time to,ti, ...,tjv-i with the time interval At. The total observation time is T = NAt and N 3> 1. 
Then, we can make N p = MiM — l)/2 pulsar pairs from M pulsars. We denote the timing residual 
of i-th pulsar at time t a as Ri(t a ). The correlation coefficient of z-th pair is defined as 



ri^K^RiMRkita), (37) 



a=0 



where %\ and %2 is the number allocated to the first and second pulsar in the z-th pair. Under the 
existence of the isotropic stochastic GW background, the ensemble average of T{ is [H7| 



(r i ) = o s q {(e i ), (31 



where a g is the root mean square of the timing residuals induced by the GW background and given 
by 

■2 



a 9 



h Pn(f)df. (39) 



The highest and lowest frequency of GWs are given by fh = 1/2 At and fi = 1/T. The power 
spectrum of the timing residuals Pr is defined as (Ri(f)Ri(f')) = Pji(f)6(f — f')/2 and that induced 
by the GW background is given by 



PrU) = g^Gw(/)- (40) 

Note that induced by the GW background has a specific dependence on the angle between the 
direction to ii-th and zVth pulsars, Q{. This dependence is characterised by ({0i), which is described 

as 

. ... 3 , x 1 _ . . 1 - cos 9 , 1 ; for x = 

((9) = -xlnx - - + -(1 + S(x)), x = ,6(x) = < . . (41) 



4 2 2 ; otherwise 
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The signature of the GW background can be extracted by calculating the following quantity, 



where f and ( are the arithmetic mean over all pairs of pulsars, and of and are the sample variance 
of r and (. We define this quantity as a "signal" . If pulsars are distributed isotropically, ( = and 
a\ — \j a/48. Hereafter, we replace ( and by these values. The ensemble average of S is given by 



2 

(S) = , ^ , (43) 



where 



*L = ^E<^-<^>) 2 >- (44) 

p i=0 

If there is no correlation between data of different pulsars, (rj) = 0, S follows a Gaussian distribution 
with zero mean and variance 1/N P . Therefore, we define the SNR as p = (S) WN P . 

We assume that the noise of each pulsar is white and uncorrelated with that of the other pulsar 
and the signal of GWs. We also assume that all pulsars have the same noise and denote the root 
mean square of time residuals induced by noise as a n . Then, the SNR becomes 



M(M — 1)/2 



1 + [x(l + <?) + 2K/<7 9 ) ! + (aja,)'] /Na 2 ' 



where 

^ N-1N-1 

* = E E c '&' c <^ = (Hta)Ri(h)) ■ (46) 

9 a=0 6=0 

Assuming that Pr(J) is the monotonically decreasing function of /, the SNR can be enhanced 
by low-pass filtering and whitening, which modifies the SNR to [89J 



; m(m -i)/2 

" = ,/ TTiL7W (47) 



where 

^max 7— » / ■ \ ^ A N, 



~ 2 2^ ^ ( P ^YcA (48) 



~ 2 

i=l " v y i=l 



Here, we define the discrete power spectra for the z-th frequency bin as 

"fi+l r.fi+1 Ofy 

dfP R (f), P n {i) = \ 
'fi Jfi 



rfi+i r.fi+i 9rr 

P g (i) = / dfP R (f), P n (i) = / dfP n (f) = P d (i) = P g (i) + P B (i), (49) 
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where 



fi 



; for i 



-0.5 



for i > 1 



(50) 



The summation is carried out only over the frequency bins in which the GW signal dominates the 
noise, P g (i) > P n , and we define N max as the number of the highest frequency bin. 

Since we use the discrete power spectra and iV max is discrete, the SNR given in Eq. (JIT]) is a 
discontinuous function of model parameters. This is inevitable as long as the TOA data are sampled 
at discrete time intervals. However, for the Fisher matrix calculation, we replace the equations by 
the expression of the continuous power spectrum, 




fn 



^2 



(J 



fh ~ fl Jfi 
C /max 



df 



PrU) 



N(f h -fi)J fl 



df 



Pn(f) 



c 2 



(51) 

(52) 
(53) 



Here, / max is given by P R (f mekX ) = P„(/ max ). Then, Eq. 
parameters, and is a good approximation of Eq. ff4T|) for iV n 
the Fisher matrix as 

1 8S 8S 



(|5ip become a smooth function of the 
lx 3> 1. Using Eq. (jBTT) . we can calculate 



N 2 86, 86, 



(54) 



where S = (l + a^Jo^af) 1/2 and N = (M(M — l)/2) If the interval of the discontinuity 
in Eq. (|4"T|) is much smaller than the error width derived from the Fisher matrix Eq. (I54p . it is 
consistent to approximate Eq. (|47p by Eq. (]5ip and Eq. f|S4]) provides a good prediction. This is 
the case in our fiducial model investigated in the next section. 

In this paper, we consider SKA, the future radio telescope array which is expected to discover 
a large number of pulsars and observe pulses with high accuracy, and take the parameters as N = 
500, M = 100, <j n = 50ns, T = lOyears, according to Ref. [88]. The detection threshold is taken as 
p > 4 again. NANOGrav, one of the latest experiments, have placed the upper bound of the GW 
background, GW < 1.9 x 10~ 8 for / ~ l/(5years) = 6.3 x 10~ 9 Hz, under the assumption that GW 
has a power-law spectrum [30] . 



-1/2 



3.4 Measurement of CMB fluctuation 

The Fisher matrix for measurement of CMB fluctuation is given by 



l X,X' 1 3 



(55) 
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bands [GHz] 


#FWHM 


arcmin 




a T [/iK] 


CTp [/iK] 


70 


14.0 


4.7 


6.7 


100 


10.0 


2.5 


4.0 


143 


7 


1 


2.2 


4.2 


217 


5 





4.8 


9.8 



Table 1: Survey parameters adopted in our analysis for Planck. The values are taken from Ref. |47j . 



bands [GHz] 


#FWHM 


arcmin 




a T [/iK] 


CTp [/iK] 


45 


17 


5.85 


8.27 


70 


11 


.0 


2.96 


4.19 


100 


8 





2.29 


3.24 


150 


5 


2.21 


3.13 


220 


3 


5 


3.39 


4.79 



Table 2: Survey parameters adopted in our analysis for CMBpol. The values are taken from Ref. [4"8] . 



where X and X' are summed over the temperature(TT), E-mode polarization(EE), B-mode polariza- 
tion(BB) and cross-correlation between temperature and E-mode(TE). Cov is the covariance matrix 
and given by 



Cov(Cf T ,Cf T ) 



Cov(C EE ,C EE ) 



cov(c?r,cn 



(2l+l)fs 



(cr T +N T ,r, 

(c* E +N P ,r, 
{cr+N P ,)\ 
[{erf + (cr + N T ,)(cr + n p ,)] 



(21 






2 


(21 






2 


(21 





TT r~<EE\ 



cov(ar,c l 



TT s-iTE\ 



(2/ + l)/ ( 



{C 



TE\2 



Cov(Cr,C< 



Cov(Cf* 



(2/+1)/, 
2 



cr E (cr+N P ,), 



n ' (2l + l)f s 
Cov(Cr,Cf B ) = Cov(Cf E ,Cf B ) = Cov(C? E ,C BB ) 



0. 



(56) 
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where f s denotes the sky coverage and is set 0.65 for both Planck and CMBpol in this paper. The 
noise power spectrum N TiP -i is given by 



T,P;i 



/ 



Nr 



T.Pd 



7 FWHM°T,P 



exp 



1(1 



(0 

FWHM 



81n2 



(57) 



where 0^p WHM is the full width at half maximum of the Gaussian beam and a^'p is the root mean 
square of the instrumental noise per pixel for temperature or polarization, for the i-th frequency 
band. The frequency bands and parameter values are provided in Tables [1] and [2] for Planck and 
CMBpol, respectively. 

The current limit on Gfi from CMB experiments is derived in Ref. [91] using WMAP and SPT 
data [90J. Their analysis provides the limit of f str < 0.0175, where 



a 



str 



str 



a 



a 



str 



inf 



2000 . 

E ZI t 1 n TT,str 



1=2 



in 



inf 



2000 . 

i- 1 riTT,inf 



i=2 



4tt 



(5* 



P 



We find this constraint corresponds in our model to G7/i < 1.4 x 10 7 for p = 1, Gfi < 3.6 x 10 8 for 
lO" 1 and Gfi < 1.0 x 10" 8 for p = 10" 2 . 
Future experiments will improve the limit for G[i by measuring the B-mode. The minimum value 
of Gfi reachable by the B-mode measurement can be estimated by yjcr^Gn) 2 [3H], where 



f.T, 



21 + 1 



BB,res 



+ N 



P\i 



BB,str 



(Gfif 



(59) 



and C, ' res is the residual noise of Cf B after removing the contamination from the foregrounds. For 



the case where the lensing effect is not removed, c^ B ' res = (jf B,len , we can probe Gfi > 2.4 x 10 8 
for p = 1, Gfi > 8.0 x 10~ 9 for p = 10 _1 , and Gfi > 2.6 x 10^ 9 for p = 10~ 2 by Planck, and 
Gfi > 1.2 x 10^ 8 for p= 1, Gfi> 3.9 x 10~ 9 for p = 10 _1 , and Gfi > 1.2 x 10~ 9 for p = 10" 2 by 
CMBpol. Even if we consider the case where the lensing effect is perfectly removed, (jf B)res = 0, the 
above values are not improved significantly, because the instrumental noise we assume here is larger 
than the lensing noise, cf B,len -C Np.j. 



iBB.len 



4 Constraints on string parameters from future experiments 

In this section, we forecast constraints on the cosmic string parameters expected from various types of 
future experiments, using the Fisher matrix calculations. We investigate three fiducial models where 
different types of experiments complement to determine the parameters with better accuracy. Here, 
we include only string parameters Gfi, a, p in theoretical parameters and calculate the Fisher matrix 

4 Note that the string network model assumed in Ref. [HI] is different from that in this paper. This causes a small 
difference in the shape of the CMB spectrum. We neglect this difference and apply their constraint on f tsr to our 
model. 
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for them. Assuming that the cosmological parameters are determined with sufficient accuracy, we set 
them to the aforementioned values and do not marginalize the likelihood over them when calculating 
constraints on the string parameters. Before that, we show the accessible parameter space by current 
and future experiments. 

4.1 Accessible parameter space of cosmic string search 



Figs. 4(a) , 4(b) , and 4(c) show the regions in the a — G\i plane which are excluded by current 
experiments or cosmological constraints and those which can be probed by future experiments for 
p — 1, 10 _1 , and 10 -2 . We also show the fiducial parameter points which we consider in the following 
Fisher analysis. 

For GW burst, the region above each curve represents the parameter space where bursts from 
cosmic strings are detectable more than once per year by each detector. The threshold amplitude 
and the most sensitive frequency /best for each interferometer are fh = 1.7 x 10 -21 and /best = 150Hz 
for current LIGO, fh = 3.4 x 10" 23 and / bcs t = 220Hz for Advanced LIGO, fh = 5.8 x 10" 21 and 
/best = 7.0 x l(T 3 Hz for eLISA/NGO, and fh = 1.9 x 1(T 24 and / best = 0.25Hz for BBO/DECIGO. 
For GW background, the detection threshold is determined whether the amplitude of the background 
flew exceeds the sensitivity of the detector at / bcst , which is GW = 7.2 x 1CT 6 for current LIGO, 
n GW = 4.5 x 1(T 9 for Advanced LIGO, and i2 G w = 9.2 x 1(T 17 for BBO/DECIGO. We assume 
that if the lowest frequency of GWs emitted by strings, fi c ~ (at ) _1 , where t is the present age of 
the universe, is lower than /best; the GWs cannot be detected either as a background or burst. The 
vertical cutoff on the left side of the curves for eLISA/NGO and BBO/DECIGO correspond to the 
value of a which gives f lc = / best . 

For the latest pulsar timing experiment, we use the constraint on flaw from NANOGrav. The 
parameter space is excluded if strings predict the GW background larger than flew = 1-9 x 1CT 8 at 
/ = l/(5years). Here, we assume that NANOGrav cannot detect GWs if f c > l/(5years), which 
corresponds to the vertical line at the left. For SKA, we show the region where SNR exceeds 4. 

We also show the cosmological constraints from CMB and BBN. The constraint is derived from 
the fact that the energy density of the GW background must be small at the last scattering and BBN, 
so as not to distort the fluctuations of the CMB or not to change abundance of various nuclei. The 
CMB constraint is J f2cw(/)^l n / < 1-4 x 10~ 5 at the last scattering [96] and the BBN constraint 
is / n GW (f)d\nf < 1.6 x 1(T 5 at the epoch of BBN [131197]. The lower limit of the integral is 
determined by the lowest frequency of the GWs emitted by largest and youngest loops at the time 
of CMB and BBN. The upper limit is the frequency of GWs emitted by the earliest loops, which we 
assume to be formed at the end of the friction domination, when the temperature of the Universe is 
~ \fG\i. 



From Figs. 4(a) , 4(b) and 4(c) , we find that, for large a, current and future pulsar timing experi- 
ments are powerful to search for cosmic strings. On the other hand, pulsar timing experiments cannot 
access to small a, and other types of experiments, such as the direct detection by interferometers 
or CMB measurement help to set constraints. In the following subsections, we investigate future 



5 If we take into account the fact that each interferometer has sensitivity over some frequency range, GWs may be 
detectable even if /';,. > /best- This may slightly expand the accessible region for large Gfx. 
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a 



(c) p = 1Q- 2 

Figure 4: Accessible parameter space in the G\i-a plane for p — 1, p — 1CT 1 and p = 1CT 2 . The 
colored regions are excluded by current experiments or cosmological considerations. The region 
above the solid or dashed lines can be probed by each future experiment. For GW direct detection 
experiments, solid lines correspond to a background search, and dashed lines correspond to burst 
detection. Here, "Adv. LIGO+" means the future interferometer network consisting of Advanced 
LIGO, Advanced Virgo and KAGRA and "current pulsar" means the limit from NANOGrav. The 
numbered crosses denote the fiducial points studied in Sees. I4.2H4.4I 
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(a) The burst rate at / = 220Hz. The bursts in the 
blue region are detectable as rare bursts by the future 
ground-based interferometer network. 
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Figure 5: The burst rate at f — 220Hz and the background spectrum Oqw for G/j, = 10 9 , a 
I0~ 9 ,p= 1. 



constraints on the cosmic string parameters from the different types of experiments by choosing 



parameter sets indicated in Fig. 4 (a 



4.2 Case 1: GfJ, = HT 9 ,a = 10~ 9 ,p = 
and pulsar timing experiments 



1 — ground-based interferometers 



First, we consider the case for G\i = 10~ 9 ,a = 10~ 9 ,p = 1. In this case, a is large enough for 
GWs to be detected in the frequency range of pulsar timing experiments. For such a value of a, 
the tension is already severely constrained and Gp = 10 -9 is the maximum value allowed by current 
pulsar timing constraints. Because of the small tension, we do not expect detection of the string 
signature by future CMB experiments. However, we instead expect future pulsar timing experiments 
such as SKA and the ground-based interferometers will detect GWs from strings. 

For this parameter set, future interferometers can detect 168 rare bursts with p > 4, where we 
assume a 3-year run of the interferometer network consisting of Advanced LIGO, Advanced Virgo 
and KAGRA. However, the GW background is not large enough to be detected at / ~ 10 2 Hz. In 
contrast, a 10-year observation by SKA can detect the GW background with p ~ 33. We show 
the burst rate dR/dlnh estimated at / = 220Hz, the most sensitive frequency of the ground-based 
interferometers, as a function of burst amplitude in Fig. 5(a) We also show the spectrum Qqw m 
Fig. 5(b) The plateau region of the spectrum seen in Fig. 5(b) corresponds to GWs emitted in the 



radiation-dominated era, and the bump in the low frequency region corresponds to those emitted 
after matter-radiation equality. Since the energy density of GWs is diluted compared with the total 
energy density after matter-radiation equality, the GWs emitted in the radiation-dominated era is 
more suppressed compared to that emitted recently, and this makes background detection at high 
frequencies difficult. 

we show the expected constraints on the string parameters from 



In Figs. 6(a), 6(b) and 6(c 



the burst detection by the ground-based interferometers and the GW background detection by SKA, 
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log a 
(c) a — p 



Figure 6: Marginalized 2a constraints on cosmic string parameters shown in the G/j — a, Gji — p 
and a — p planes. The fiducial parameter set, denoted by the black cross, is taken to be Gfi = 10~ 9 , 
a = 10~ 9 , p = 1. The solid black line represents the constraints from the burst detection alone by the 
ground-based interferometer network consisting of Advanced LIGO, Advanced Virgo and KAGRA. 
The red line represents the combined constraints from the burst detection and the GW background 
measurement by SKA. 
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estimated by Fisher matrix calculations. Since the constraints from SKA alone are quite weak, we 
only show the combined constraints in the figure. One can see SKA slightly improve the constraints 
when it is combined to the constraint from the ground-based interferometers. This is because the 
constraints from pulsar timing and interferometer experiments have different directions of parameter 
degeneracy. 

Let us briefly discuss the parameter degeneracy. In this parameter set, a < TGp, so loops 
evaporate soon after its formation by emitting GWs of frequency / ~ (at)" 1 , where t is the time of 
the GW emission. We provide the rough estimates of the parameter degeneracy in dR/ dlnh and Oqw 
for a < TGp in Appendix A. For burst detection, the bursts whose amplitude is comparable to the 
sensitivity of the ground-based interferometers, fh ~ 3.4 x 10~ 23 , is in the range of < h < h 3>3 of 
Eq. (IBTJl) . where fh 3 ^ — 2.3 x 10 -27 and fh 33 = 2.2 x 10 -22 in this case. So the parameter degeneracy 
of the burst rate is oc (G/i) 3//8 a~ 3//4 p _1 . Here, h 32 and /i 3i3 are the corresponding amplitude to 
characterize when the bursts are emitted (For details, see Appendix A or Ref. |20j). For example, 
bursts who has amplitude of h 3t 2 < h < h 3>3 are emitted between the matter-radiation equality and 
z ~ 1. The rough estimate of the background spectrum is given by Eq. (IBTj) . For the GW background 
measurement by SKA, the second term dominates in Eq. (!6Tl) . so the parameter degeneracy is 



In this fiducial model, eLISA/NGO can also detect GW bursts from strings. Fig. 7(a) shows 
the burst rate at the most sensitive frequency, / = 7 x 10~ 3 Hz. In the case of eLISA/NGO, the 
detectable bursts corresponds to the case of h > h 3t3 of Eq. fl60|) . where fh 3t3 = 6.9 x 10~ 21 for 
/ = 7 x 10 _3 Hz. Such bursts are emitted recently at z < 1. In Figs. 7(b), 7(c), and 7(d), we show 



the constraints expected from a 3-year run of eLISA/NGO, which makes 1.4 x 10 4 burst detections 
with p > 4. One can clearly see eLISA/NGO can provide much stronger constraints than ground- 
based interferometers and SKA. And, of course, BBO/DECIGO will determine the parameters with 
a significant accuracy. Our Fisher calculation indicates the expected errors on the parameters are 
0(0.1)%. 

4.3 Case 2: Gfi = 10 ~ 7 , a = 10 16 , j> =1 — ground-based interferometers 
and CMB experiments 

Next, we study the case where Gfi = 10~ 7 , a = 10~ w ,p = 1. In this case, a is extremely small, 
so pulsar timing experiments cannot detect GWs from strings. This means that the tension is not 
constrained strongly by the current pulsar timing experiments and, if G[i ~ 10~ 7 which is still 
allowed by current CMB and LIGO experiments, we can expect future CMB experiments to find 
string signatures. Also, we can expect both burst and background detection by future ground-based 
interferometers. The ground-based interferometers will detect 1.8 x 10 5 rare bursts with p > 4 and the 
GW background with p ~ 187, where we again assume a 3-year run of the world-wide interferometer 



network. We show the burst rate at / = 220Hz in Fig. 8(a) and the background spectrum f2cw in 



Fig 8(b) 



In Figs. 6(a), 6(b) and 6(c), we show the constraints on the string parameters from the CMB 



6 Note that, the direction of the degeneracy seen in the figures does not directly corresponds to the parameter 
dependence described here, since the shown constraints are marginalized over the other parameter. 
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log G\i log a 

(c) G/i-p (d) a— p 

Figure 7: (a): The burst rate at / = 7 x l(r 3 Hz for G\i = 10~ 9 ,a = 10~ 9 ,p = 1. (b),(c),(d): 
Marginalized 2a constraints from eLISA on cosmic string parameters in G/i — a, Gfi — p and a — p 
planes, respectively. The black cross represents the fiducial point. 
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(a) The burst rate at / = 220Hz. The bursts in the 
blue region are detectable as rare bursts by the future 
ground-based interferometer network. The bursts in 
the orange region consists a GW background. 
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Figure 8: Same as Fig. |5]but for Gfi = 10' 
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P = 1. 



observation and the burst and GW background detection by the ground-based interferometers. For 
the GW background, we use the weak signal approximation, Eq. (1361) . to calculate the Fisher matrix. 
For CMB, we consider the constraints from Planck and CMBpol. We derive the constraints from 
CMBpol neglecting the lensing effect. The results are not significantly affected by including the 
lensing effect, as discussed in Sec. 3.4. 

In this case, loops are extremely short-lived. The bursts whose amplitude is comparable to the 
sensitivity of the ground-based interferometers, fh ~ 3.4 x 10 -23 , corresponds to h > /i 33 in Eq. (IPO]) , 
where fh 33 = 4.7 x 10~ 25 in this case. These bursts are emitted recently at z < 1. Therefore, for 



-l 



the burst detection, the direction of the parameter degeneracy is oc (Gp) a ' p . The background 



spectrum is again expressed by Eq. (IfjTl) . The bump-like spectrum in Fig. 8(b) corresponds to the 



second term in Eq. (I6T1) . which represents GWs emitted recently, and its parameter dependence is 
oc Gpa~ 1 ^ 3 p~ 1 . CMB measurements provide information on infinite strings, which is characterized 
by only Gfi and p, and do not contain information on a. The overall amplitude of CMB spectra is 
proportional to (G/i) 2 and decreasing p leads to enhancement of the amplitude and the change of the 
spectral shape as explained in Sec. 3. We numerically find that the dominant contribution to the 
Fisher matrix comes from the temperature spectrum around 1000 < / < 2000. For such values of /, 
QTT.str - g j-Qygjjjy proportional to p~ 2 , as shown in Fig. 
of the CMB constraint is oc (Gjj) 2 p~ 2 . 



2(b) Therefore, the parameter degeneracy 



4.4 Case 3: G\i = 10 
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Finally, we study the case where Gfi = 10" , a = 10 ,p = 1. For this small value of G/i, only 
BBO/DECIGO can detect string signals. In this case, BBO/DECIGO detects 35 rare bursts with 
p > 4 in a 3-year run and measure the GW background with very high SNR, p ~ 510. We show 
the burst rate estimated at the best frequency of BBO/DECIGO, / = 0.25Hz, in Fig. 10(a) and the 
background spectrum in Fig. 10(b) 



23 




-16.3 -16.2 -16.1 -16 -15.9 -15.8 -15.7 

log a 

(c) a — p 



Figure 9: Marginalized 2a constraints on cosmic string parameters shown in the Gjj — a, Gji — p 
and a — p planes. The fiducial parameter set, denoted by the black cross, is taken to be G/j, = 10 -7 , 
a = 10~ 16 , p — 1. The solid black line represents the constraints from the burst detection by the 
interferometer network consisting of Advanced LIGO, Advanced Virgo and KAGRA alone. The red 
line represents the combined constraints from the burst detection and the GW background measure- 
ment by the interferometer network. The blue or green lines are the constraints from combination 
of interferometers and Planck or CMBpol, respectively. The orange lines are the constraint from 
CMBpol alone. 
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(a) The burst rate at / = 0.25Hz. The bursts in the 
blue region are detectable as rare bursts by the future 
ground-based interferometer network. The bursts in 
the orange region consists a GW background. 

Figure 10: Same as Fig. [5] but for G\i = 10 



10" 4 10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 



frequency[Hz] 

(b) The background spectrum f^cw- 



-11 



a 



10 



-13 



p — 1. 



Figures 11 (a 



11(b) and 11(c) show the expected constraints on the string parameters from the 
burst detection and the background measurement by a 3- year run of BBO/DECIGO. Here, we use 
the exact formula of the Fisher matrix for the GW background, Eq. ( 135]) , because of the high SNR. 
The background measurement provides stronger constraints than burst detection. However, the 
constraints from background measurement intrinsically have strong parameter degeneracies, since 
the information is practically only one flow & t / = /best- Thus, although constraints from the burst 
detection is weak because of the small number of detectable events, it dramatically tightens the 
errors when combined with the constraints from the background detection. This is again thanks to 
the difference in the parametric dependencies. 

In this parameter set, loops are marginally short-lived. The bursts detectable by BBO/DECIGO 
are emitted recently, z < 1, and their rate is expressed by the case of h > h^^ of Eq. flBTJj) . where 
/^3,3 — 4.5 x 10~ 28 . So the parameter degeneracy is oc (G/x) 2 a 1/,3 p _1 . The background spectrum 
around / = 0.25Hz corresponds to GWs emitted at z < 1 and is roughly expressed by the second 



term of Eq. 061|) . whose parameter dependence is oc Gpia 



-1/3 



V 



5 Summary 

Among many types of cosmic string signatures, gravitational waves from cosmic string loops and 
CMB fluctuations induced by infinite strings are important and future experiments will help to test 
the existence of cosmic strings. If detected, we may even able to extract information on the nature 
of cosmic strings. In this paper, we extend our previous work which investigated the constraints on 
the string parameters from future ground-based GW experiments. In addition to the ground-based 
interferometers, this paper have investigated constraints from space-borne interferometers, pulsar 
timing arrays and CMB experiments. Furthermore, we have studied the combination of information 
from these observations. 
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Figure 11: The 2a constraints from BBO/DECIGO on cosmic string parameters in the G/i — a, 
Gfi — p and a — p planes. The black line shows the marginalized 2a constraints from background 
measurement alone. The red line represents the combined constraints from burst detection and 
background measurement. The fiducial parameter set, denoted by the black cross, is taken to be 
Gn = 10" 14 , a = 10" 13 , p= 1. 
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Each experiment sheds light on the different aspect of cosmic strings and gives us different in- 
formation on strings. More technically speaking, the constraints from different experiments have 
different parameter degeneracies. CMB experiments probe infinite strings through their gravita- 
tional effects on the background photons, while GW interferometers and pulsar timing experiments 
probe string loops by detecting GWs emitted from them. Among GW experiments, each of them 
targets a different frequency band; ground-based interferometers detect GWs of / ~ 220Hz, space- 
borne ones probe those of / ~ 10~ 4 to 0.1Hz and pulsar timing experiments are sensitive to those of 
/ ~ 10~ 8 Hz. GWs of different frequency are emitted at different epoch of the Universe and provide 
us with independent information. Besides, different types of GW observations, the burst detection 
and the background measurement, provide different information on cosmic strings. These are the 
reasons why we can break the parameter degeneracies by combining these experiments and obtain 
better constraints on cosmic string parameters. 

In this paper, we have studied three different fiducial models, where different types of experiments 
are help each other to constrain the string parameters. The first case is G[i = 10 _9 ,a = 10~ 9 ,p = 
1, where both future GW interferometers and pulsar timing such as SKA can detect GWs from 
cosmic strings. We also calculate constraints from eLISA/NGO and found that it is more efficient 
to constrain string parameters than the pulsar timing and ground-based experiments. The second 
case is Gfi = 10~ 7 , a = 10~ 16 ,p = 1, where both future GW interferometers and CMB experiments 
can detect string signatures. The third case is G/i = 10~ 14 , a = 10~ 13 ,p = 1, where the tension is so 
small that only ultimate space-borne interferometers such as BBO and DECIGO can detect string 
signatures. We have shown that future GW interferometers, especially space-borne ones, are very 
powerful to investigate cosmic strings not only because of their extreme sensitivity but also because 
they can probe strings in two different ways, the background measurement and the burst detection. 

Finally, we should note that, in the case where we can determine the parameters with a very 
good accuracy, uncertainties in the string network model become more important. In that case, 
further theoretical study will be needed to perform analysis with more accurate modeling, or precise 
measurements of cosmic string GWs by future space-borne interferometers may even be able to shed 
light on the theoretical uncertainties in the string evolution. 

Appendix 

A Parameter dependence of the burst rate and the GW 
background spectrum 

In this appendix, we provide rough estimates of the burst rate dR/d\nh and the GW background 
spectrum f2 G w to understand the dependence of these quantities on the cosmic string parameters. 
We concentrate on the case of a < TGfi, where loops are short-lived, since this case applies to all 
the fiducial models investigated in Sec. 4. For other cases, see the appendices of Ref. [20] . 
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In this case, a < TG/i, the burst rate is given by 

' (GW^a-VV 2 (fy) " n/1 ° r l8 /% 12/5 h^ 

(G/x) 3 / 8 a- 3 / 4 7 - 2 /- 5/2 io 17/8 ^- 11/8 " (60) 

;/l 3 ,2 < h < G/Xa 2 /3/-4/3 t -l/3 = ^ 
(G / x) 2 a 1 /3 7 -2 / -14/3 t -8/3^- 3 . h > ^ 

Bursts in the range of hz \ < h < h 3 2 , h 3 2 < h < /13 3, and fa > ^3 are emitted in the radiation- 
dominated era, in the matter-dominated era, and at z <§C 1, respectively. The burst rate for h < h 3jl 
is suppressed strongly 

The background spectrum is given by 

n GW (/) ~ g^;^ + Gry-ta-wtffif-w (61) 

Here, the first term is the contribution from GWs emitted in the radiation-dominated era and the 
second one represents GWs emitted recently. Note that the former has smaller amplitude than the 
latter by the factor of Q r /Q m , since the energy density of GWs decays faster than the total energy 
density in the matter-dominated era. 

Acknowledgment 

K.M. and T.S. would like to thank the Japan Society for the Promotion of Science for financial 
support. This work is supported by the Grant-in-Aid for Scientific Research from the Ministry of 
Education, Science, Sports, and Culture (MEXT), Japan, No. 23340058 (S.K.), No. 24740149 (S.K.), 
No. 23740179 (K.T.), No. 24111710 (K.T.) and No. 24340048 (K.T.). 

References 

[1] T. W. B. Kibble, J. Phys. A 9, 1387 (1976). 

[2] A. Vilenkin and E. P. S. Shellard, "Cosmic Strings and Other Topological Defects," Cambridge 
University Press, Cambridge, England (1994). 

[3] S. Sarangi and S. H. H. Tye, Phys. Lett. B 536, 185 (2002) |arXiv:hep-th/0204074] . 

[4] N. T. Jones, H. Stoica, S. H. H. Tye, Phys. Lett. B563, 6-14 (2003). [hep-th/0303269) . 

[5] G. Dvali and A. Vilenkin, JCAP 0403, 010 (2004) [arXiv:hep-th/0312007] . 

[6] A. Vilenkin, Phys. Lett. B107, 47-50 (1981). 

[7] C. J. Hogan, M. J. Rees, Nature 311, 109-113 (1984). 



OR 
din h 



28 



R. R. Caldwell and B. Allen, Phys. Rev. D 45, 3447 (1992). 

[9] R. R. Caldwell, R. A . Battye and E. P. S. Shellard, Phys. Rev. D 54, 7146 (1996) 
[arXiv:astro-ph/9607T30| . 

[10] T. Damour and A. Vilenkin, Phys. Rev. Lett. 85, 3761 (2000) [arXiv:gr-qc/0004075] . 

[11] T. Damour and A. Vilenkin, Phys. Rev. D 64, 064008 (2001) |arXiv:gr-qc/0104026| . 



[12] T. Damour and A. Vilenkin, Phys. Rev. D 71, 063510 (2005) |arXiv:hep-th704 10222|. 

[13] X. Siemens, V. Man dic and J. Creighton, Phys. Rev. Lett. 98, 111101 (2007) 
|arXiv:astro-ph/0610920] . 

[14] M. R. DePies and C. J. Hogan, Phys. Rev. D 75, 125006 (2007) |arXiv: astro-ph/0702335| ; 

[15] M. Kawasaki, K. Miyamoto, K. Nakayama, Phys. Rev. D81, 103523 (2010). [ arXiv:1002.0652l 
[astro-ph.CO]]. 

[16] S. Olmez, V. Mandic, X. Siemens, Phys. Rev. D81, 104028 (2010). [ arXiv:1004.0"890l [astro- 
ph.CO]]. 

[17] P. Binetruy, A. Bohe, C. Caprini and J. -F. Dufaux. larXiv:1201.0983l [gr-qc]. 

[18] P. Binetruy, A. Bohe, T. Hertog and D. A. Steer, Phys. Rev. D 82, 126007 (2010) 
|arXiv:1009.2484l [hep-th]]. 

[19] A. Bohe, Phys. Rev. D 84, 065016 (2011) jarXiv: 1103.07681 [hep-th]]. 

[20] S. Kuroyanagi, K. Miyamoto, T. Sekiguchi, K. Takahashi and J. Silk, arXiv: 1202.3032] [astro- 
ph.CO]. 

[21] D. Sigg [LIGO Scientific Collaboration], Class. Quant. Grav. 25, 114041 (2008). 

[22] T. Accadia et al, Class. Quant. Grav. 28, 114002 (2011). 

[23] K. Kuroda [LCGT Collaboration], Class. Quant. Grav. 27, 084004 (2010). 

[24] P. Amaro-Seoane, S. Aoudia, S. Babak, P. Binetruy, E. Berti, A. Bohe, C. Caprini and M. Colpi 
et al, larXiv:1201.362T1 [astro-ph.CO]. 

[25] E. S. Phinney et al. 2003, Big Bang Observer Mission Concept Study (NASA). 

[26] S. Kawamura, M. Ando, N. Seto, S. Sato, T. Nakamura, K. Tsubono, N. Kanda and T. Tanaka 
et al, Class. Quant. Grav. 28, 094011 (2011). 

[27] R. N. Manchester, AIP Conf. Proc. 983, 584 (2008) [ arXiv:0710.5"026l [astro-ph]]. 

[28] R. D. Ferdman, R. van Haasteren, C. G. Bassa, M. Burgay, I. Cognard, A. Corongiu, N. D'Amico 
and G. Desvignes et al, Class. Quant. Grav. 27, 084014 (2010) |arXiv:1003.3lQ5l [astro-ph. HE]]. 



29 



F. Jenet, L. S. Finn, J. Lazio, A. Lommen, M. McLaughlin, I. Stairs, D. Stinebring and J. Ver- 
biest et al. larXiv:0909.1058l [astro-ph.IM]. 

P. B. Demorest, R. D. Ferdman, M. E. Gonzalez, D. Nice, S. Ransom, I. H. Stairs, Z. Arzouma- 
nian and A. Brazier et al. larXiv7l 201.6641 [astro-ph.CO]. 



http://www. skatelescope.org 



N. Kaiser and A. Stebbins, Nature 310, 391 (1984). 
J. R. Gott, III, Astrophys. J. 288, 422 (1985). 

U. Seljak, U. -L. Pen, N. Turok, Phys. Rev. Lett. 79, 1615-1618 (1997). |astro-ph/970423l] . 
L. Pogosian and T. Vachaspati, Phys. Rev. D 60, 083504 (1999) [astro-ph/990336T] . 
K. Benabed, F. Bernardeau, Phys. Rev. D61, 123510 (2000). 

L. Pogosian, S. H. H. Tye, I. Wasserman and M. Wyman, Phys. Rev. D 68, 023506 (2003) 
[Erratum-ibid. D 73, 089904 (2006)] [ hep-th/0304188] . 

U. Seljak, A. Slosar, Phys. Rev. D74, 063523 (2006). |astro-ph/060ll~43] . 

N. Bevis, M. Hind marsh, M. Kunz and J. Urrestilla, Phys. Rev. Lett. 100, 021301 (2008) 
|astro-ph/070'2223| [ASTRO-PH]]. 



N. Bevis, M. Hindmarsh, M. Kunz, J. Urrestilla, Phys. Rev. D76, 043005 (2007). 
|arXiv:0704.3"800 [astro-ph]]. 

L. Pogosian, M. Wyman, Phys. Rev. D77, 083509 (2008). |arXiv:0711.0"747l [astro-ph]]. 

M. Kawasaki, K. Miyamoto and K. Nakayama, Phys. Rev. D 82, 103504 (2010) |arXiv:1003.370T1 
[astro-ph.CO]]. 

R. Battye and A. Moss, Phys. Rev. D 82, 023521 (2010) [ arXiv:1005.0479l [astro-ph.CO]]. 

N. Bevis, M. Hindmarsh, M. Kunz, J. Urrestilla, Phys. Rev. D82, 065004 (2010). 
[arXiv: 1005.26631 [astro-ph.CO]]. 

D. Yamauchi, K. Takahashi, Y. Sendouda, C. -M. Yoo and M. Sasaki, Phys. Rev. D 82, 063518 
(2010) [arXiv:1006.0687l [astro-ph.CO]]. 

P. Mukherjee, J. Urrestilla, M. Kunz, A. R. Liddle, N. Bevis, M. Hindmarsh, Phys. Rev. D83, 
043003 (2011). [arXiv: 1010.56621 [astro-ph.CO]]. 



[Planck Collaboration], astro-ph/0604069 



D. Baumann et al. [CMBPol Study Team Collaboration], AIP Conf. Proc. 1141, 10 (2009) 
[arXiv:0811.3"9T9l [astro-ph]]. 



30 



[49] D. P. Bennett, F. R. Bouchet, Phys. Rev. Lett. 60, 257 (1988), Phys. Rev. Lett. 63, 2776 (1989), 
Phys. Rev. D41, 2408 (1990). 



[50 
[51 

[52 
[53 
[54 
[55 
[56 
[57 

[58; 

[59 
[60 

[61 

[62; 

[63 
[64 

[65; 

[66 
[67 



B. Allen, E. P. S. Shellard, Phys. Rev. Lett. 64, 119-122 (1990). 

G. R. Vincent, M. Hindmarsh, M. Sakellariadou, Phys. Rev. D56, 637-646 (1997). 
|astro-ph/9612135] . 

V. Vanchurin, K. D. Olum, A. Vilenkin, Phys. Rev. D74, 063527 (2006). |gr-qc/0511159] . 

C. Ringeval, M. Sakellariadou, F. Bouchet, JCAP 0702, 023 (2007). |astro-ph/05lT646] . 
C. J. A. P. Martins, E. P. S. Shellard, Phys. Rev. D73, 043515 (2006). [astro-ph/0511~792] . 
K. D. Olum, V. Vanchurin, Phys. Rev. D75, 063521 (2007). |astro-ph/06104l9] . 

X. Siemens, K. D. Olum, A. Vilenkin, Phys. Rev. D66, 043501 (2002). [gr-qc/0203006] . 



J. Polchinski, J. V. Rocha, Phys. Rev. D74, 083504 (2006). [hep-ph/0606205|. 

F. Dubath, J. Polchinski, J. V. Rocha, Phys. Rev. D77, 123528 (2008). [ arXiv:0711.0"994l [astro- 
ph]]. 

V. Vanchurin, Phys. Rev. D77, 063532 (2008). |arXiv:0712.2"236 [gr-qc]]. 

L. Lorenz, C. Ringeval and M. Sakellariadou, JCAP 1010, 003 (2010) [arX iv:1006.093T1 [astro- 
ph.CO]]. 



http: //www.sfu.ca/~levon/cmbact.html 



S. Foreman, A. Moss and D. Scott, Phys. Rev. D 84, 043522 (2011) [arXiv: 1106.40181 [astro- 
ph.CO]]. 

A. Avgoustidis, E. P. S. Shellard, Phys. Rev. D73, 041301 (2006). | astro-ph/05 12582] . 

K. Takahashi, A. Naruko, Y. Sendouda, D. Yamauchi, C. -M. Yoo, M. Sasaki, JCAP 0910, 003 
(2009). |arXiv:0811.4698l [astro-ph]]. 

C. J. A. P. Marti ns, E. P. S. Shellard, Phys. Rev. D54, 2535-2556 (1996). 
[arXiv:hep-ph/9602271| . 

C. J. A. P. Martins and E. P. S. Shellard, Phys. Rev. D 65, 043514 (2002) [hep-ph/0003'298"| . 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) [arXiv: 1001 .45381 
[astro-ph. CO]]. 



[68] X. Siemens, J. Creighton, I. Maor, S. Ray Majumder, K. Cannon and J. Read, Phys. Rev. D 
73, 105001 (2006) |gr-qc/0603li5] . 



31 



[69] B. P. Abbott et al. [LIGO Scientific Collaboration], Phys. Rev. D 80, 062002 (2009) 
|arXiv:0904.47T8l [astro-ph.CO]]. 



70 
71 

72 
73 

74 
75 

76 
77 
78 

79 

80 
81 

82 
83 
81 
85 
86 

87 



M. Tegmark, A. Taylor and A. Heavens, Astrophys. J. 480, 22 (1997) |arXiv:astro-ph/960302l] . 

M. I. Cohen, C. Cutler and M. Vallisneri, Class. Quant. Grav. 27, 185012 (2010) jarXiv: 1002.41531 
[gr-qc]]. 

J. S. Key and N. J. Cornish, Phys. Rev. D 79, 043014 (2009) jarXiv: 08 12 .1590 [gr-qc]]. 

N. J. C ornish and L. J . Rubbo, Phys. Rev. D 67, 022001 (2003) [Erratum-ibid. D 67, 029905 
(2003)] |gr-qc/02090iT] . 

B. Abbott et al. [LIGO Collaboration], Phys. Rev. D 69, 102001 (2004) |gr-qc/0312056] . 

L. Baggio et al. [AURIGA and LIGO Scientific Collaborations], Class. Quant. Grav. 25, 095004 
(2008) |arXiv:0710.0497l [gr-qc]]. 

L. S. Finn and A. N. Lommen, Astrophys. J. 718, 1400 (2010) jarXiv: 1004 .34991 [astro-ph.IM]]. 

B. Allen and J. D. Romano, Phys. Rev. D 59, 102001 (1999) |arXiv:gr-qc/9710117] . 

H. Kudoh, A. Taruya, T. Hiramatsu and Y. Himemoto, Phys. Rev. D 73, 064006 (2006) 
|gr-qc/051H45] . 

A. Nishizawa, A. Taruya, K. Hayama, S. Kawamura and M. a. Sakagami, Phys. Rev. D 79, 
082002 (2009) jarXiv: 0903 .05281 [astro-ph.CO]]. 

V. Corbin and N. J. Cornish, Class. Quant. Grav. 23, 2435 (2006) |arXiv:gr-qc/05 12039] . 

B. Abbott et al. [LIGO Scientific Collaboration], Phys. Rev. D 69, 122004 (2004) 
|gr-qc/03120~88] . 

N. Seto, Phys. Rev. D 73, 063001 (2006) |arXiv:gr-qc/05 10067] . 

F. B. Estabrook and H. D. Wahlquist, Gen. Relativ. Graviation, 6, 439 (1975) 

M. V. Sazhin, Soviet Astron., 22, 36 (1978) 

S. Detweiler, Astrophys. J. 234, 1100 (1979) 

D. R. Lorimer and M. Kramer, "Handbook of Pulsar Astronomy," Cambridge University Press, 
Cambridge, England (2005). 

R. W. Hellings and G. S. Downs, Astrophys. J. 265, L39 (1983) 

A. Sesana and A. Vecchio, Class. Quant. Grav. 27, 084016 (2010) [ arXiv:1001.3T6Tl [astro- 
ph.CO]]. 



32 



[89] F. A. Jenet, G. B. Hobbs, K. J. Lee and R. N. Manchester, Astrophys. J. 625, L123 (2005) 
|astro-ph/0504458| . 

[90] R. Keisler, C. L. Reichardt, K. A. Aird, B. A. Benson, L. E. Bleem, J. E. Carlstrom, C. L. Chang 
and H. M. Cho et al, Astrophys. J. 743, 28 (2011) jarXiv: 1 1 05 . 31821 [astr o-ph. CO] ] . 

[91] C. Dvorkin, M. Wyman and W. Hu, Phys. Rev. D 84, 123519 (2011) [arXiv:1109.4947l [astro- 
ph.CO]]. 

[92] A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 538, 473 (2000) lastro-ph/99m77] . 



[93] http: / / camb.info/ 



[94] M. Zaldarriaga, D. N. Spergel and U. Seljak, Astrophys. J. 488, 1 (1997) |[ astro-ph/97021~57| . 
[95] L. Knox, Phys. Rev. D 52, 4307 (1995) [| astro-ph/9504054| . 

[96] T. L. Smith, E. Pierpaoli, M. Kamionkowski, Phys. Rev. Lett. 97, 021301 (2006). 
[astro-ph/0603144] . 

[97] R. H. Cyburt, B P. Fields, K. A. Olive, E. Skillman, Astropart. Phys. 23, 313-323 (2005). 
[astro-ph/0408033] . 



33 



